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^ (57) Abstract: The present invention relates to processing of spread spectrum signals, where a continuous signal of a comparatively 
© high frequency is received (1610). This signal is sampled at a basic sampling rate whereby a resulting sequence of time discrete 
® signal samples is produced, which are in turn quantised into a corresponding level-discrete sample value (1620). A plurality of data 

words are formed, which each includes one or more consecutive sample values (1630). Information obtained from these data words 
^ is correlated (1640) with at least one representation of a signal source specific code sequence, which has been pre-generated in the 

form of a code vector (1600). The correlation step specifically involves correlating at least each vector in a sub-group of the code 
^ vectors with at least one vector that has been derived from the data word. Thereby resulting data is produced (1650). 
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Spread Spectrum Signal Processing 

5 THE BACKGROUND OF THE INVENTION AND PRIOR ART 

The present invention relates generally to processing of spread 
spectrum signals by means of vector based algorithms. More 
particularly the invention relates to a method of receiving spread 
spectrum signals according to the preamble of claim 1 and a 
10 signal receiver according to the preamble of claim 28. The 
invention also relates to a computer program according to claim 
26 and a computer readable medium according to claim 27. 

Spread spectrum technology is becoming increasingly important 
in many areas, of which cellular communication systems and 

15 global navigation satellite systems (GNSS) represent two 
important examples. Moreover, the large variety of transmission 
standards creates a demand for hybrid- or general-purpose 
receivers. For instance, both the cellular standards cdma2000 
and WCDMA involve transmission of spread spectrum signals. 

20 However, due to differences in the air interfaces one and the 
same terminal cannot be used in the two systems. Instead, a 
dedicated terminal must be employed in each system. Alter- 
natively, a dual mode terminal must be designed, which includes 
two separate transceiver chains. Correspondingly, a navigation 

25 receiver adapted for one GNSS, say the Global Positioning 
System (GPS; U.S. Government), is not able to receive signals 
from a satellite that belongs to a different GNSS, such as the 
Galileo system (the European programme for global navigation 
services) or the Global Orbiting Navigation Satellite System 

30 (GLONASS; Russian Federation Ministry of Defense). 
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In order to make this possible, a multi-mode receiver must be 
designed. However, including multiple receiver chains in a 
single device is not only expensive, it also renders the unit bulky 
and heavy, particularly if more than two signal formats are to be 
5 processed. Therefore, a programmable software receiver 
solution is desired whose signal processing principles may be 
altered according to which signals that shall be received and 
processed. 

Various software solutions are already known for processing 
10 GNSS signals. The patent document WO02/50561 is one 
example in which a GPS receiver tracking system for guiding 
missiles is described. The receiving station here includes one or 
more processing stages for receiving and processing GPS data. 
A control signal controls whether a quantity of processing stages 
15 should be increased or decreased in response to the result of a 
correlation operation. The tracking and recovery of timing in- 
formation is entirely performed in software. However, the pro- 
cessing involves real-time calculation of relatively complex 
Fourier transforms, and therefore requires considerable pro- 
20 cessing resources. Thus, a receiver of this type is not parti- 
cularly suitable for small sized units, such cellular phones or 
portable GNSS-receivers. 

The article "TUPM 12.4: Software Solution of GPS Baseband 
Processing", ULSI Laboratory, Mitsubishi Electric Corporation, 

25 IEEE, 1998 by Asai T., et al. outlines the principles for a soft- 
ware GPS receiver using an embedded microprocessor, which 
enables signals from up to eight satellites to be demodulated at 
around 40 MIPS (millions of / or Mega Instructions Per Second). 
A 2 MHz GPS baseband signal is here fed to the software radio 

30 receiver, which performs spectrum despreading and mixing; 
synchronous demodulation, satellite acquisition and tracking; as 
well as phase compensation. An algorithm is proposed through 
which the intermediate frequency (IF) signal is 1-bit sampled 
and where downconversion is performed before despreading. 

35 Resulting data in the form of 32 samples are continuously 
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processed in parallel. A DRAM (Dynamic Random Access 
Memory) stores a particular C/A (Coarse Acquisition) code table 
for each satellite from which signals are received. All other data 
and instructions are stored in a cache memory in order to 
5 accomplish a high-speed processing. Although this strategy is 
declared to result in a reasonably efficient implementation, the 
document lacks a specific description as to how the dispreading 
and the correlation should be effected in order to, in fact, obtain 
this effect. 

10 Akos. D. et al., "Tuning In to GPS - Real-Time Software Radio 
Architectures for GPS Receivers", GPS World, July 2001 
describes a receiver architecture through which IF signal 
samples are fed directly from a radio front-end to a program- 
mable processor for continued processing. The article mentions 

15 the possibility of using single instruction multiple data (SIMD) 
instructions to process multiple data samples in parallel. How- 
ever, portability- and flexibility problems are recognized, and 
again, there is no teaching as to how an efficient parallel pro- 
cessing could actually be accomplished. 

20 Dovis, F. et al., "Design and Test-Bed Implementation of a 
Reconfigurable Receiver for Navigation Applications", Electro- 
nics Department, Politecnico di Torino, Navigation Signal Ana- 
lysis and Simulation Group, Spring of 2002 relates to the design 
of a reconfigurable GNSS receiver which is capable of fusing 

25 data from two or more different GNSS:s. The document sketches 
an architecture which, in addition to a radio front-end, includes a 
Field Programmable Gate Array (FPGA) and a Digital Signal 
Processor (DSP). The authors address various computational 
load issues. However, they do not present an algorithm which 

30 fulfills the identified requirements. 

Hence, the prior art includes a number of examples of software- 
based GNSS-receivers. Nevertheless, there is yet no distinct 
teaching of a highly efficient solution which is suitable for imple- 
mentation in software and has capabilities that are at least in par 
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with those of today's ASIC-based solutions for receiving and pro- 
cessing spread spectrum signals in real-time (ASIC = Application 
Specific Integrated Circuit). 

SUMMARY OF THE INVENTION 

5 The object of the present invention is therefore to provide a 
solution for receiving and processing spread spectrum signals, 
which solves the problems above and thus to present a truly 
software adapted strategy that may be implemented in a 
physically small and power efficient device. 

10 According to one aspect of the invention the object is achieved 
by a method of receiving spread spectrum signals as initially 
described, which is characterized by comprising a preparation 
for the correlation step, wherein the preparation takes place 
before the receiving the continuous signal. This preparation 

15 involves pre-generating a multitude of code vectors, which each 
represents a particular code sequence of the at least one signal 
source specific code sequence. Moreover, according to the 
invention, the correlation step Involves multiplying at least each 
vector in a sub-group of the code vectors with at least one 

20 vector, which is derived from the data word. 

This strategy is advantageous because the proposed vector 
approach makes it possible for a digital processor (e.g. a micro- 
processor) to process multiple signal samples in parallel during 
each clock cycle and thus make very efficient use of the 
25 processor. Furthermore, by pre-generating the code vectors 
valuable processing capacity is saved, which in turn renders it 
possible to deal with a sufficient amount of input data per time 
unit in order to, for example track navigation signals that are 
transmitted by the satellites in a navigation satellite system. 

30 According to a preferred embodiment of this aspect of the 
invention, each code vector represents a particular signal 
source specific code sequence, which is sampled at the basic 
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sampling rate and quantised with the quantising process that is 
used to produce the level-discrete sample values. Preferably, 
the code vectors are also further adapted to the format of the 
incoming data signal in a way which renders the following 
5 processing efficient. 

According to another preferred embodiment of this aspect of the 
invention, the at least one signal source specific code sequence 
represents so-called pseudo random noise. This type of code 
sequences is beneficial in environments where two or more 

10 different signals are modulated onto the same carrier frequency. 
Namely, the different pseudo random noise signals are ortho- 
gonal (or at least almost orthogonal) to each other. Conse- 
quently, a first information signal spread by means of a first 
pseudo random noise signal causes a minimal deterioration of & 

15 second information signal which is spread by means of a second 
pseudo random noise signal, and vice versa. 

According to another preferred embodiment of this aspect of the 
invention, the receiving step involves down conversion of an 
incoming high-frequency signal to an intermediate frequency 

20 signal. The high-frequency signal is presumed to have a spec- 
trum which is symmetric around a first frequency and the inter- 
mediate frequency signal is presumed to have a spectrum which 
is symmetric around a second frequency, which is considerably 
lower than the first frequency. Hence, the receiving step 

25 transforms the information carrying signal to a baseband which 
is technically more easy to handle than the band at which the 
original incoming signal is located. 

According to another preferred embodiment of this aspect of the 
invention, the method includes the following steps. First, a 
30 maximum frequency variation of the second frequency due to 
Doppler-effects is determined. Then, a Doppler frequency interval 
around the second frequency is defined. The Doppler frequency 
interval has a lowest frequency limit equal to the difference 
between the second frequency and the maximum frequency 
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variation, and a highest frequency limit equal to the sum of the 
second frequency and the maximum frequency variation. After that, 
the Doppler frequency interval is divided into an integer number of 
equidistant frequency steps, and subsequently a frequency 
5 candidate vector is defined for each frequency step. This way of 
dividing the spectrum is advantageous because it provides a 
flexible modeling of the intermediate carrier frequency and its 
variations. 

According to another preferred embodiment of this aspect of the 
10 invention, an integer number of initial phase positions is . 
determined for the frequency candidate vector. The integer 
number thus represents a phase resolution, such that a 
relatively high number corresponds to a comparatively high 
phase resolution, whereas a relatively low number corresponds 
15 to a comparatively low phase resolution. Furthermore, a carrier 
frequency-phase candidate vector is defined for each combi- 
nation of carrier frequency candidate vector and initial phase 
position. The different frequency-phase candidate vectors are 
beneficial, since they allow the determination of a very accurate 
20 estimate of frequency- and phase characteristics of a received 
signal, which in turn vouches for a demodulated signal having a 
high quality. 

According to yet another preferred embodiment of this aspect of 
the invention, the number of elements in each carrier frequency- 

25 phase candidate vector is determined. Subsequently, the carrier, 
frequency-phase candidate vectors are stored according to a 
data format, which is adapted to performing multiplication 
operations between the data word and a segment of a carrier 
frequency-phase candidate vector. This adaptation is 

30 advantageous because thereby several signal samples may be 
processed in parallel by means of comparatively simple 
operators. Preferably, the adapting of the data format involves 
adding at least one element to each segment of the carrier 
frequency-phase candidate vector, such that the segment 

35 attains a number of elements which is equal to the number of 
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elements in each of the at least one vector that is derived from 
the data word. Hence, it is namely possible to process a 
segment of the carrier frequency-phase candidate vector to- 
gether with one of the at least one vector by either a so-called 
5 Single Instruction Multiple Data (SIMD)-operation or a logical 
XOR-operation. 

According to still another preferred embodiment of this aspect of 
the invention, the method also includes the following steps. 
First, a maximum variation of the code rate due to Doppler- 

10 effects is determined. Then, a Doppler rate interval is defined 
around a center code rate. The Doppler frequency interval has a 
lowest code rate limit equal to the difference between the center 
code rate and the maximum code rate variation, and a highest 
frequency limit equal to the sum of the center code rate and the 

15 maximum code rate variation. In analogy with the intermediate 
frequency spectrum, the Doppler rate interval is divided into an 
integer number of equidistant code rate steps, and a code rate 
candidate is defined for each code rate step. In further analogy 
with the intermediate frequency, an integer number of possible 

20 initial code phase positions is determined for each code rate 
candidate. The integer number thus represents a code phase 
resolution, such that a relatively high number corresponds to a 
comparatively high code phase resolution, whereas a relatively 
low number corresponds to a comparatively low code phase 

25 resolution. Typically, the spread spectrum code modulation is 
less sensitive to distortion caused by Doppler effects than the 
frequency modulation. Therefore, the code rate steps may be 
relatively large. However, a certain degree of Doppler shift esti- 
mation is desirable to attain a demodulated signal of high quality. 

30 According to yet another preferred embodiment of this aspect of 
the invention, a set of code vectors is generated for each signal 
source specific code sequence by sampling each code rate- 
phase candidate vector with the basic sampling rate, and thus 
produce a corresponding code vector. Again, this simplifies a 

35 parallel processing of several signal samples by means of 
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comparatively simple operators. 

According to still another preferred embodiment of this aspect of 
the invention, a modified code vector is generated on basis of 
each code vector by: copying a particular number of elements 
5 from the end of an original code vector to the beginning of the 
modified code vector, and copying the particular number of 
elements from the beginning of the original code vector to the 
end of the code vector. This extension of the code vector is 
highly advantageous because it allows extraction of tracking 
10 parameters through early-, prompt- and late techniques by 
performing a simple translation of local copies of the signal 
source specific code sequence. 

According to still another preferred embodiment of this aspect of 
the invention, a set of modified code vectors is stored for each 

15 signal source specific code sequence. Each modified code 
vector here contains a number of elements, which represent a 
sampled version of at least one full code sequence. As 
mentioned above, the modified code vector also includes a 
repetition of the beginning and the end of the sequence. A 

20 particular modified code vector is defined for each combination 
of code rate candidate and code phase position. These pre- 
generated modified code vectors save considerable processing 
capacity in the processor, since they need not be calculated in 
real time, but may simply be acquired from a memory. 

25 According to yet another preferred embodiment of this aspect of 
the invention, the data format of the modified code vectors is 
adapted with respect to the data format of the at least one 
vector that has been derived from the data word, such that a 
modified code vector and one of the at least one vector may be 

30 processed jointly by a SIMD-operation or an XOR-operation. 
Naturally, this is advantageous from a processing efficiency 
point-of-view. 

According to yet another preferred embodiment of this aspect of 
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the invention, the method involves an initial acquisition phase 
and a subsequent tracking phase. The acquisition phase esta- 
blishes a set of preliminary parameters which are required for 
initiating a decoding of signals that are received during the 
5 tracking phase. The parameters may include: a modified code 
vector, a carrier frequency candidate vector an initial phase 
position, a code phase position and code index, which denotes a 
starting sample value for the modified code vector. A successful 
acquisition phase thus results in that at least one signal source 
10 specific code sequence is identified and that this signal is 
possible to track by the receiver (i.e. may be continued to be 
received). 

According to still another preferred embodiment of this aspect of 
the invention, the tracking in turn, involves the following steps. 

15 First, based the tracking characteristics, a prompt pointer is 
calculated for each modified code vector. The prompt pointer 
indicates a start position for the code sequence. The initial 
prompt pointer value is set equal to the code index. Then, at 
least one pair of early- and late pointers is assigned around 

20 each prompt pointer. The early pointer specifies a sample value 
which is positioned at least one element before the prompt 
pointer's position, and correspondingly, the late pointer specifies 
a sample value being positioned at least one element after the 
prompt pointer's position. These pointers are used to maintain 

25 the tracking of the received signal by a repeated repositioning of 
the pointers around a correlation maximum value, such that the 
prompt pointer is positioned as close as possible to this value, 
the early pointer is positioned in time somewhat prior to this 
value and the late pointer is positioned in time somewhat after 

30 this value. An important advantage attained by this strategy is 
that one and the same vector may be re-used in order to 
represent different delays. Thus, valuable memory space and 
processing capacity is saved. 

According to still another preferred embodiment of this aspect of 
35 the invention, the tracking involves the following further steps. 
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First, a sequence of incoming level-discrete sample values is 
received. Then, data words are formed from the sample values, 
such that each data word contains a number of elements which 
is equal to the number of elements in each carrier frequency- 
5 phase candidate vector. After that, a relevant set of carrier 
frequency-phase candidate vectors is calculated for the data 
word, and a pre-generated in-phase representation respective a 
quadrature-phase representation of the vector is acquired for 
each carrier frequency-phase candidate vector in the relevant 

10 set. Subsequently, each data word is on one hand multiplied 
with a in-phase representation of the carrier frequency-phase 
candidate vector in the relevant set to produce a first inter- 
mediate-frequency-reduced information word, and on the other 
hand multiplied with a quadrature-phase representation of the 

15 carrier frequency-phase candidate vector in the relevant set to 
produce a second intermediate-frequency-reduced information 
word. The intermediate-frequency-reduced information "words 
thus constitute a respective component of the vectors mentioned 
above that have been derived from the incoming data words. 

20 This multiplied strategy is advantageous because it enables a 
very efficient usage of the capacity of a general-purpose 
microprocessor. 

According to yet another preferred embodiment of this aspect of 
the invention, XOR- or SIMD-operations are used to perform the 
25 multiplication between the data word and the in-phase repre- 
sentation of the carrier frequency-phase candidate vector re- 
spective between the data word and the a quadrature-phase 
representation of the carrier frequency-phase candidate vector. 

According to another preferred embodiment of this aspect of the 
30 invention, the tracking involves the further steps of: correlating 
the first intermediate-frequency-reduced information word with a 
modified code vector starting at a position indicated by the 
prompt pointer to produce a first prompt-despread symbol string, 
correlating the first intermediate-frequency-reduced information 
35 word with a modified code vector starting at a position indicated 
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by the early pointer to produce a first early-despread symbol 
string, correlating the first intermediate-frequency-reduced infor- 
mation word with a modified code vector starting at a position 
indicated by the late pointer to produce a first late-despread 
5 symbol string, correlating the second intermediate-frequency- 
reduced Information word with a modified code vector starting at 
a position indicated by the prompt pointer to produce a second 
prompt-despread symbol string, correlating the second inter- 
mediate-frequency-reduced information word with a modified 

10 code vector starting at a position indicated by the early pointer 
to produce a second early-despread symbol string, and corre- 
lating the second intermediate-frequency-reduced information 
word with a modified code vector starting at a position indicated 
by the late pointer to produce a second late-despread symbol 

15 string. Preferably, a resulting data word is then derived for each 
set of the despread symbol strings. This may, for example be 
performed either by performing the relevant adding operation, or 
by looking up a pre-generated value in a table based on the bit 
pattern of the despread symbol strings depending on whether 

20 the resulting data words contain packed or un-packed 
information. Data of high quality may thus be obtained with a 
minimum of real time calculations, which is advantageous from a 
processing point-of-view. 

According to yet another preferred embodiment of this aspect of 
25 the invention, either XOR- or SIMD-operations are used also to 
perform the multiplication between the intermediate-frequency- 
reduced information words and the modified code vectors, since 
again, this is advantageous from a processing point-of-view. 

According to another preferred embodiment of this aspect of the 
30 invention, certain pieces of information are propagated in 
connection with completing the processing of a current data 
word and initiating the processing of a subsequent data word. 
Preferably, this propagated information includes: a pointer which 
indicates a first sample value of the following data word, a group 
35 of parameters that describe the relevant set of carrier 
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frequency-phase candidate vectors, the relevant set of code 
vectors, and the prompt-, early-, and late pointers. Hence, based 
on the propagated information, the subsequent data word may 
be processed immediately. 

5 According to a further aspect of the invention the object is 
achieved by a computer program directly loadable into the 
internal memory of a computer, comprising software for perfor- 
ming the above proposed method when said program is run on a 
computer. 

10 According to another aspect of the invention the object is 
achieved by a computer readable medium, having a program 
recorded thereon, where the program is to make a computer 
perform the above proposed method. 

According to another aspect of the invention the object is 
15 achieved by the signal receiver as initially described, which is 
characterized in that the proposed computer program is loaded 
into the memory means, such that the receiver will operate 
according to the above-described method. 

Consequently, the invention. offers,..an excellent instrument for 
20 receiving and decoding any type of spread spectrum signals in a 
general software radio receiver. For instance, one and the same 
receiver may be used to track navigation satellite signals of 
different formats. These signals may even be received in 
parallel, such that a first signal is received on a first format 
25 while one or more other signals are received on a second or 
third format. 

Moreover, the proposed solution permits one and the same 
receiver to be utilized for completely different purposes, such as 
communication in a cellular system. 

30 BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention is now to be explained more closely by 
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means of preferred embodiments, which are disclosed as 
examples, and with reference to the attached drawings. 

Figure 1 shows a spectrum graph over an exemplary spread 
spectrum signal from which information may be ex- 
tracted according to the invention, 

Figure 2 shows a spectrum graph over a frequency down- 
converted signal corresponding to the spectrum 
shown in figure 1, 

Figure 3 is an enlarged version of figure 2 in which a 
Doppler shifting of the spectrum is illustrated, 

Figure 4 illustrates how a data signal is modulated onto a 
signal source specific code sequence according to 
an embodiment of the invention, 

Figure 5 depicts a three-dimensional graph over carrier 
frequency-phase candidate vectors to be used as 
initial correlators according to an embodiment of 
the invention, 

Figure 6 shows an alternative representation the carrier 
frequency-phase candidate vectors in figure 5, 

Figures 7-8 show graphs which illustrate how a code index and 
a code phase are defined for a received signal 
according to an embodiment of the invention, 

Figures 9a, b illustrate how a modified code vector is generated 
based on an original code vector according to an 
embodiment of the invention, 



Figure 10 



30 Figure 11 



depicts a three-dimensional graph over code rate- 
phase candidate vectors to be used as subsequent 
correlators according to an embodiment of the 
invention, 

illustrates how a set of code sequence start 
pointers is assigned to a modified code vector 
according to an embodiment of the invention, 
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Figure 12 illustrates how data words are formed from a 
sequence of incoming level-discrete sample values 
according to an embodiment of the invention, 

Figure 13 illustrates how the data words of figure 12 are 
5 correlated with various pre-generated vectors 

according to an embodiment of the invention, 

Figure 14 illustrates, how multiple pairs of pointers may be 
assigned according to an embodiment of the 
invention, 

10 Figure 15 shows a proposed signal receiver, and 

Figure 16 illustrates, by means of a flow diagram, the general 
method of processing spread spectrum signals 
according to the invention. 



DESCRIPTION OF PREFERRED EMBODIMENTS OF THE 
15 INVENTION 

In order to realize a truly software based signal processing 
solution in real-time at high sampling rates, such as those 
required in a GNSS, the architecture of the microprocessor 
system must be utilized very efficiently. The present invention 

20 uses an approach, which involves a form of parallel processing 
based on pre-generated code vector representations. The 
invention thereby strikes a balance between processing speed 
and memory usage, which compensates for the capacity deficit 
of a general purpose microprocessor in relation to an ASIC. The 

25 details of this solution will become apparent from the below 
disclosure. 

Figure 1 shows a spectrum graph over an exemplary spread 
spectrum signal from which information may be extracted 
according to the invention. For example, the signal may be a 
30 GPS C/A-signal in the L1-band with a carrier frequency f HF of 
1,57542 GHz. This means that the signal's spectrum is 
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symmetrical around the carrier frequency f H F- Moreover, the 
signal has at least 98% of its energy distributed within a 2 MHz 
wide frequency band BW H f around the carrier frequency f H F- 
Therefore, it is generally sufficient if the spread spectrum signal 
5 is downconverted to an intermediate frequency f| F at around 1 
MHz in order to enable further signal processing at a lower 
frequency. 

Figure 2 shows a spectrum graph over such a frequency down- 
converted signal, whose spectrum is symmetrical around an 

10 intermediate frequency f| F , which is considerably lower than the 
carrier frequency f hf> for example f| F = 1,25 MHz. Due to move- 
ments of the transmitter (i.e. typically a satellite in orbit) and the 
receiver (i.e. typically a mobile GNSS-receiver) a Doppler shift 
may occur in the received signal, which is proportional to the 

15 relative velocity between the transmitter and the receiver. A 
maximum Doppler shift f D is defined which in the GNSS case 
normally is approximately 10 kHz. Thus, the entire spectrum 
may be shifted within an interval 2f D , such that the center 
frequency falls any where between fiF+fp (corresponding to a 

20 maximum positive Doppler shift) and f| F -fo (corresponding to a 
maximum negative Doppler shift). For illustrative purposes the 
interval 2f 0 , has been severely exaggerated in the figure 2. The 
figure also indicates a phase variation interval 90 (along an 
independent axis q>) around the intermediate frequency f iF . The 

25 phase variation interval <p D demonstrates the fact that the phase 
position of the received signal is initially unknown to the 
receiver, and thus constitutes a calibration parameter. The 
effects of this parameter will be further elucidated below with 
reference to figure 5. 

30 An enlarged version of the spectrum in figure 2 is shown in 
figure 3. The intermediate spectrum may, due to possible 
Doppler effects, be shifted such that its center frequency lies in 
the interval 2f D from f| F . m i n =fiF-fD to fiF-max=fiF + fo. The potential 
maximum shifted spectra are here indicated with a respective 

35 dashed line. The diagram in figure 3 also shows a Nyquist 
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frequency r s /2 corresponding to a minimal sampling rate r s 
which, when used for sampling the intermediate frequency 
signal, will result in a non-aliased discrete spectrum. Provided 
that f 0 = 10 kHz, f, F = 1,25 MHz, BW HF = 2 MHz and the inter- 
5 mediate frequency signal is lowpass filtered through an analog 
filter having a 2,5 MHz wide passband, the minimal sampling 
rate r s becomes 5 MHz (i.e. two times the highest frequency 
component after the analog filter). 

Figure 4 illustrates schematically how a data signal D is 

10 modulated onto a signal source specific code sequence CS on 
the transmitter side according to an embodiment of the inven- 
tion. According to a preferred embodiment of the invention, the 
signal source specific code sequence CS constitutes so-called 
pseudo random noise. The data signal D here contains a data 

15 symbol sequence [+1, -1, -1, +1, -1], and has a relatively low 
symbol rate, say 50 Hz. The signal source specific code 
sequence CS, however, has a relatively high symbol rate (or 
more correctly chipping rate). For instance, a signal source 
specific code sequence CS in the form of a GPS C/A code may 

20 have a chipping rate of 1,023 MHz and contain 1023 chips per 
period. Each chip, ch, is either +1 or -1. Hence, the C/A code 
repeats itself one per ms. The data signal D is modulated onto 
(or spread by) the signal source specific code sequence CS by 
multiplying each data symbol with the code sequence CS. The 

25 data symbol +1 thereby results in an unaltered code sequence 
CS, whereas the data symbol -1 results in an inverted code 
sequence CS. 

If a data signal D having a rate of 50 Hz is spread by means of a 
signal source specific code sequence CS having a chipping rate 
30 of 1,023 MHz, this results in 20 entire code sequences CS per 
data symbol. Namely, the period time for one data symbol is 20 
ms, whereas the period time for the code sequence is only 1 ms. 

Figure 5 depicts a three-dimensional graph over different carrier 
frequency-phase candidate vectors. As mentioned above, the 
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intermediate frequency spectrum may be Doppler shifted. 
Therefore, the particular frequency of the received signal is 
initially unknown within a Doppler frequency interval from fiF-min to 
fiF-max- Moreover, the phase position of the received signal is 
5 beforehand unknown to the receiver. In order to enable an 
estimation of the actual intermediate carrier frequency, the 
Doppler frequency interval fiF-min - fiF-max is divided into an integer 
number of equidistant frequency steps. Of course, the larger the 
number of steps is, the more accurately the intermediate carrier 

10 frequency may be estimated. In this example, a step size Af of 
20 Hz is chosen. In some applications, however, a much larger 
step size Af may provide a sufficient frequency accuracy, whereas 
in other applications (where relatively high accuracy is needed) a 
considerably smaller step size Af may be required. In any case, a 

15 Doppler frequency interval f| F -min - fiF-max of 20 kHz and a step size 
Af of 20 Hz requires 1000 steps, where a first step includes the 
frequencies from fiF-min to f| F . m in + Af, a second step includes the 
frequencies from f| F . m in + Af to fiF-min + 2Af, and so on up to a last step 
including the frequencies from f| F . m in + 999Af to fiF-max- 

20 A frequency candidate vector f| F -c is defined for each frequency 
step. In order to estimate an initial phase position of the 
intermediate carrier frequency, an integer number of initial phase 

positions 9 0 <P 7 are determined. The illustrated example 

shows eight such positions <P C . This gives a phase resolution of 

25 rc/4. For each frequency candidate vector f ( F-c. eight different 
phase positions <P 0 <P 7 are thus defined which each repre- 
sents a particular phase shift from 0 to 7n/4. A combination of a 
certain carrier frequency candidate vector f|F-c and a particular 
initial phase position <P C is referred to as a carrier frequency- 

30 phase candidate vector. Figure 5 shows all such vectors (i.e. 
eight) for the first step, where f| F -c = fiF-mm. and the second step, 
where f, F -c = fiF-mm+Af. 

Additionally, figure 5 includes an axis L c along which the 
number of sample values (or length) of the frequency candidate 
35 vector f| F -c is indicated. The L c number is preferably a power of 2 
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because this generally renders the implementation more simple. 
The number of sample values in each frequency candidate vector 
fi F . c is set in consideration of the available memory space, the 
desired frequency resolution and the desired phase resolution. If, 
5 in the example above, a length of 512 samples is chosen and each 
sample requires one byte, the total (uncompressed) memory 
demand would be 3,9 MB (i.e. 1000 frequencies x 8 phase 
positions x 512 B = 4096000 B = 3,9 MB). 

Figure 6 shows an alternative representation the carrier 
10 frequency-phase candidate vectors V f(p (fiF-c. ^c) in figure 5. 
Here, the vectors V f(p (f, F .c, <P C ) range from the carrier frequency 
fiF-c=fiF-min (say 1,24 MHz) to f| F -c=fiF-max (say 1,26 MHz), from 
the phase shift position <P c =( Po (say 0) to <P C =( P7 (say 7n/4), all 
have a length L c of s n samples (say 512) and are represented in 
15 a cube. Hence, a particular frequency-phase candidate vector 
Vf<p(fiF.c. ^c) »s given by a cuboid segment in the cube which is 
parallel with the L c -axis, and ranges from L c =Si to L c =s n . 

After having generated all the carrier frequency-phase candidate 
vectors V f)p (fiF-c» ^c). the vectors are stored in a digital memory 

20 with a comparatively short access time, such that they may be 
readily accessed during a later acquisition and/or tracking 
phase. Preferably, the carrier frequency-phase candidate 
vectors V f()) (fiF.c. ^c) are stored according to a data format which 
is adapted to performing multiplication operations between 

25 incoming data words (representing one or more received 
signals) and segments of carrier frequency-phase candidate 
vectors V f(p (f| F .c, <Pc)- The vectors may either include packed 
samples or they may contain samples that are unpacked to a 
smallest manageable unit in the microprocessor, such as a byte. 

30 Further details regarding the data format issues will be 
discussed below. 

Figure 7 shows a graph over an exemplary code sequence CS, 
which has been received, downconverted and sampled 
according to an embodiment of the invention. After subsequent 
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removal of the intermediate frequency component, the code 
sequence CS either attains the signal value +1 or -1, and the 
symbols are altered at a chipping rate of for example 1,023 
MHz. At the receiver side, the code sequence CS is sampled 

5 with a sampling rate r s of for example 5 MHz. The sampling in- 
stances are vindicated by means of dots on the code sequence 
signal CS along an axis St. As is apparent from the figure there 
are approximately four sample values per chip period. Hence, 
the requirements of the Nyquist theorem are clearly fulfilled. 

10 However, the samples are not capable of determining the 
transition instances with a sufficient accuracy, i.e. where the 
code sequence signal CS changes from representing one chip 
symbol to representing another one chip symbol. Therefore, a 
so-called code phase must be determined. Moreover, a parti- 

15 cular sample value at which the code sequence starts must also 
be established. 

Figure 8 shows a graph, which illustrates the sample value at 
which the code sequence CS starts (the so-called code index 
CI) and how a code phase, Cph, is defined in relation to this 
20 sample value. An initial segment of the code sequence signal 
CS in figure 7 is therefore represented in figure 8. 

A first sample instance is here presumed to occur shortly prior 
to the actual start of the code sequence CS signal period. 
Hence, a corresponding sample value s-i belongs to the end of 

25 the preceding period. This is illustrated by means of the line 
representing the code sequence signal CS being dashed. A 
subsequent sample instance, however, overlaps the code 
sequence CS period. The corresponding sample value s 2 is 
defined as the code index CI. The actual beginning of a code 

30 sequence signal is determined in a preceding acquisition phase 
and typically involves correlation between the received signal 
and a local copy of the code sequence CS. 

The distance between the start of the code sequence signal CS 
and the code index CI, in turn, is defined as the code phase 
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Cph, and thus represents a measure of the translation (or skew) 
between the sample instances Sj and the chip transitions. The 
sampling interval between two consecutive sample instances Sj 
is divided into an integer number of possible initial code phase 
5 positions in order to estimate the code phase Cph. In this 
example, ten 0.1 -steps ranging from 0.0 to 0.9 are defined, 
where a first step 0.0 would indicate that code sequence signal 
CS starts at the code index CI, and a last step 0.9 would 
indicate that code sequence signal CS starts almost at the 
10 preceding sample value s^ The figure 8 shows a code phase 
Cph of 0.4, i.e. a situation where the code sequence signal CS 
starts approximately in the middle between two consecutive 
sample instances s,. 

In analogy with the carrier frequency-phase candidate vectors 
15 V f< p(f| F -c. ^c). code rate-phase candidate vectors may now be 
defined for each signal source specific code sequence CS and 
combination of code rate CR, code index CI and code phase 
Cph. Nevertheless, according to a preferred embodiment of the 
invention, a modified code vector is first produced based on 
20 each original code vector. 

Figures 9a and 9b illustrate how this may be performed. A code 
vector CV in figure 9a contains a number of sample values, say 
5000. Preferably, the number of sample values is selected in 
relation to sampling rate, such that a full code sequence period 
25 is recorded. In the above example, with a sampling rate r s of 5 
MHz, this means that the 5000 samples cover exactly one full 
code sequence CS period (which has a duration of 1 ms and 
contains 1023 chips, since the chipping rate is 1,023 MHz). 

A modified code vector CV m is now generated on basis of the 
30 original code vector CV by copying a particular number of 
elements E e (e.g. two) from the end of the original code vector 
CV to the beginning of the modified code vector CV m . 
Correspondingly, a particular number of elements E b (e.g. two) 
from the beginning of the original code vector CV are copied to 
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the end of the code vector CV m> i.e. after the elements E e . 
Figure 9b shows the resulting modified code vector CV m . Since, 
in this example, each chip symbol is sampled at approximately 
four instances, two sample values are equivalent to about one 
5 half chip period. Thereby, the above copying of the elements E e 
and E b allows a one-half-chip shifting of the code vector CV 
during a correlation operation between the code vector CV and 
the received data. This will be discussed further below with 
reference to figure 11. 

10 Figure 10 depicts a three-dimensional graph over different code 
rate-phase candidate vectors CV m (CR C -c. Cph) which represent 
a signal source specific code sequence CS(i). Analogous to the 
intermediate carrier frequency, the nominal code rate CR C for 
the code rate candidate vector CR C -c may also vary due to 

15 Doppler effects from a lowest value CR min to a highest value 
CRmax- However, since the code rate (or chipping rate) is com- 
paratively high, the rate variation here becomes rather small. 

For example, a Doppler shift of 10 kHz at the intermediate 
carrier frequency is equivalent to a shift in a nominal code rate 

20 CR C at 1,023 MHz of merely 6,5 Hz. Namely, the Doppler shift in 
respect of the code rate is generally equal to the ratio between 
the carrier frequency and the code rate (1575,42 MHz / 1,023 
MHz = 1540; a Doppler shift of 10 kHz at the intermediate 
carrier frequency thus corresponds to 10000/1540 * 6,5 Hz). 

25 The code rate candidate vector CR C -c may attain a value any- 
where in a Doppler rate interval between CR min = 1022993,5 Hz 
and CR max = 102306,5 Hz. The Doppler rate interval CR min - 
CR max is divided into an integer number of equidistant code rate 
steps ACR, say 13. In this example, each interval CR mi n, 

30 CR m in+ACR CR m ax thus spans 1 Hz. According to a preferred 

embodiment of the invention, the code rate step size is adaptive, 
e.g. determined on basis of the sampling frequency and the 
desired resolution. 

In figure 10, the three-dimensional graph (or cube) contains a set 
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of code rate-phase candidate vectors CV m (CR C -c. Cph), where the 
code rate ranges from CR C -c=CR min to CR C -c=CR ma x. the code 
phase ranges from Cph=0.0 to Cph=0.9, and each vector has a 
length L C vm of s m samples (for example 5004). Hence, a particular 
5 code rate-phase candidate vector CV m (CR C -c. Cph) is given by a 
cuboid segment CV m (CR min +2ACR, 0.4) which is parallel with the 
L-cvm -axis and ranges from L C vm=Si to L C vm =s m - 

After having generated all the code rate-phase candidate 
vectors CV m (CR C -c. Cph), these vectors are stored in a digital 

10 memory with a comparatively short access time, such that they 
may be readily accessed during the acquisition and/or tracking 
phase. Preferably, the code rate-phase candidate vectors 
CV m (CRc-c. Cph) are stored according to a data format which is 
adapted to correlation operations to be performed between a 

15 vector that has been derived from incoming data words and a 
segment of a the code rate-phase candidate vectors CV m (CR C -c. 
Cph). For example, the code rate-phase candidate vectors 
CV m (CR C - C . Cph) may include packed samples or they may 
contain samples that are unpacked to a smallest manageable 

20 unit in the microprocessor, such as a byte. Further details 
regarding these data format issues will be discussed below with 
reference to figure 13. 

During tracking of a spread spectrum signal, i.e. a continued 
reception of the signal and demodulation of the data contained 

25 therein, it is necessary to update relevant tracking parameters 
for the signal. Normally, the tracking phase is preceded by a so- 
called acquisition phase during which a set of preliminary 
parameters are established that are required for initiating the 
signal decoding. A successful acquisition phase thus results in 

30 that at least one signal source specific code sequence is 
identified and a data signal transmitted by means of this 
sequence is possible to demodulate. The tracking characteris- 
tics that may be associated with a signal source specific code 
sequence include: a modified code vector CV m (which defines a 

35 carrier frequency candidate vector f| F . c and an initial phase 
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position <P C ), a code phase position, Cph, and a code index, CI. 
In order to maintain an adequate timing of the correlating 
operations in the receiver, a set of pointers are required which 
indicate the start position of the code sequence in relation to the 
5 modified code vector and are updated repeatedly, preferably 
between each correlation. 

Figure 11 illustrates one such set of pointers that are assigned 
to a modified code vector CV m . A prompt pointer P P indicates a 
current estimate of the code sequence start position. This 

10 pointer P P is initially set equal to the code index CI. Adequate 
Pp-pointer positions for any subsequent segments of the 
received data are obtained by retrieving appropriate modified 
code vectors CV m from the digital memory where such pre- 
generated are stored. Between each segment, the code rate 

15 candidate vector CR C -c as well as the initial code phase Cph is 
updated. Additionally, at least one pair of early- and late 
pointers P E and P L respectively is assigned on each side of the 
prompt pointer P P , where the early pointer P E specifies a sample 
value being positioned at least one element before the prompt 

20 pointer's P P position, and the late pointer P L specifies a sample 
value being positioned at least one element after the prompt 
pointer's P P position. Each correlation is then performed with the 
presumption that the code sequence starts at the P P -position, at 
the P E -position as well as at the P L -position. If the correlation 

25 over the P P -position results in a higher correlation value than 
over any of the other positions, and at the same time the P E - and 
P L -pointers are balanced (i.e. the correlation over these 
positions result in equal values), this is interpreted as an optimal 
setting of the tracking parameters. Should, however, this not be 

30 true the sampled data must be repositioned in relation to the P P -, 
P E - and P L -pointers, such that the correlation peak again 
coincides with the P P -pointer's position. An un-calibrated pointer 
positioning may thus be detected via a difference between the 
correlation value over the P E -pointer position and the correlation 

35 value over the P L -pointer position, while the correlation over the 
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Pp-position still results in the highest correlation value. In this 
case, the potentially highest correlation value (corresponding to 
the optimal P P -positioning) lies somewhere between the current 
Pp-position and the current P E -position or between the current 
5 P P -position and the current P L -position depending on which of 
the P E - and P L -pointer positions that results in the highest 
correlation value. A pointer updating normally takes place 
between each segment. A second or third pair of P E - and PL- 
pointers outside the first pair of pointers further enhances the 

10 possibilities of accomplishing an accurate adjustment of the 
tracking. Depending on the value of the received symbol, an 
optimal setting P P -pointer may, in fact, also be equivalent to the 
correlation over the P P -position resulting in a lower correlation 
value than over any of the other pointer positions. Namely, this 

15 is true if a negative valued data bit is received. However 
otherwise, the same principles apply. 

In any case, the P P -, P E - and P L -pointers may only be positioned 
within certain intervals in the modified code vector CV m . The P P - 
pointer may be set anywhere within a first interval 111, the P E - 

20 pointer within a second interval 112 and the P L -pointer within a 
third interval 113. All the pointers are always set together, such 
that their relative distances remain unchanged. A limit value P|_-min 
for the P L -pointer typically coincides with the first element of the 
modified code vector CV m . Namely, outside this range a 

25 meaningful correlation cannot be performed. As a consequence of 
this and the fact that the pointers are set jointly P P -, P E - and P L , a 
corresponding limit value P E . m in for the P E -pointer is also defined. 

Figure 12 illustrates how data words d(1), d(2), .... d(N) are 
formed from a sequence 1210 of incoming level-discrete sample 

30 values according to an embodiment of the invention. The data 

words d(1) d(N) are formed such that each data word d(k) 

contains a number of elements which is equal to the number of 
elements s n in each carrier frequency-phase candidate vector. For 
instance, the length of these vectors may be 512 sample values 

35 (as in the example described above with reference to the figures 
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5 and 6). Given that each sequence 1210 contains 5000 sample 
values, the number N of data words becomes 10 (5000/51 2»9, 76 
means that the last quarter of the tenth data word d(10) is empty). 
According to the invention, the elements in a current data word 

5 d(k) are processed jointly (i.e. in parallel). Then, the processing of 
a subsequent data word d(k+1) is initiated, and any processing 
parameters obtained during the processing of the current data 
word d(k) that are required when processing the subsequent data 
word d(k+1) are propagated as input data to the processing of the 

10 latter. These processing parameters preferably include: a pointer 
p d which indicates a first sample value of the subsequent data 
word d(k+1), a group of parameters which describe the relevant 
set of carrier frequency-phase candidate vectors V f(p (f| F -c. ^c) (i.e. 
a carrier frequency candidate vector f| F -c representing an in-phase 

15 version and a carrier frequency candidate vector f| F . c representing 
a quadrature-phase version), the relevant set of code vectors 
CV m , and prompt-, early-, and late pointers P P , P E , and P L 
respectively. 

Figure 13 illustrates how the data words of figure 12 are corre- 
20 lated with various pre-generated vectors according to an embo- 
diment of the invention. First, a relevant set of carrier frequency- 
phase candidate vectors V f< p(f| F - C , ^c) is calculated for the data 
word d(k). This set is normally obtained in a preceding 
acquisition phase (or, if the data word d(k) is not the first word 
25 to be processed, during the processing of a previous word). 

Thereafter, for each carrier frequency-phase candidate vector 
V f(p (fiF-c. *Pc) in the relevant set, a pre-generated in-phase repre- 
sentation f| P i and a quadrature-phase representation f (FQ of the 
vector (V f<p (f| F .c. <Pc)) respectively is acquired. Both these repre- 
30 sentations are found in the digital memory where the pre-gene- 
rated carrier frequency-phase candidate vectors V f<p (f| F .c, ^c) are 
stored because they are merely one quarter period delays of 
one and the same vector, and the memory holds initial phase 
position representations for an entire period, see figure 6. 
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Then, in order to eliminate the influence of the carrier frequency 
component, the data word d(k) is multiplied with the in-phase 
representation f ifi of the carrier frequency-phase candidate 
vector Vf 9 (f|F. c , <Pc) in the relevant set. As a result thereof, a first 
5 intermediate-frequency-reduced information word S| F . t (k) is pro- 
duced. Correspondingly, the data word d(k) is also multiplied 
with the quadrature-phase representation f| FQ of the carrier 
frequency-phase candidate vector V f<p (f| F .c, ^c) in the relevant 
set, and a second intermediate-frequency-reduced information 
10 word S| F . Q (k) is produced. 

Preferably, the multiplication between the data word d(k) and 
the in-phase representation f| F i respective the quadrature-phase 
representation f| FQ of the carrier frequency-phase candidate 
vector V f(p (f| F .c, <Pc) are performed by means of XOR- or SIMD- 
15 operations (i.e. 1-bit and multiple-bits multiplications respec- 
tively). This is namely possible if the data word d(k) and the 
vector V f(( >(f| F .c, <Pc) have compatible data formats. 

Subsequently, the first intermediate-frequency-reduced informa- 
tion word S| F .|(k) is multiplied with a modified code vector CV m . P , 

20 which has been retrieved from the digital memory where these 
pre-generated vectors are stored, and starts at a position 
indicated by the prompt pointer P P . A first prompt-despread 
symbol string A| P (k) is produced as a result of this operation. 
The first intermediate-frequency-reduced information word 

25 S, F .|(k) is also multiplied with a modified code vector CV m . E (k)» 
which starts at a position indicated by the early pointer P E . A 
first early-despread symbol string A (E (k) is produced as a result 
of this operation. The first intermediate-frequency-reduced 
information word S| F .|(k) is likewise multiplied with a modified 

30 code vector CV m . L (k), which starts at a position indicated by the 
late pointer P L , and a first late-despread symbol string An(k) is 
produced. Moreover, the second intermediate-frequency- 
reduced information word S| F . Q (k) is multiplied with a modified 
code vector CV m . P (k), which starts at a position indicated by the 

35 prompt pointer P P , and a second prompt-despread symbol string 
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Aop(k) is produced. Correspondingly, the second intermediate- 
frequency-reduced information word S| F -ci(k) is multiplied with a 
modified code vector CV m . E (k), which starts at a position 
indicated by the early pointer P E , and a second early-despread 
5 symbol string A QE (k) is produced. Likewise, the second 
intermediate-frequency-reduced information word Sip.Q(k) is 
multiplied with a modified code vector CV m . L (k), which starts at a 
position indicated by the late pointer P L , and a second late- 
despread symbol string A QL (k) is produced. 

10 After that, a resulting data word D R .| E (k), D R ., P (k), D R .| L (k), 
D R . Qe (k), D R . QP (k) and D R . QL (k) may be derived from the each of 
the despread symbol strings Aip(k), A| E (k), A| L (k), A QP (k), A QE (k) 
and A QL (k) respectively by adding the elements in the respective 
string together. If the strings represent un-packed data, the 

15 processor may simply perform the relevant adding operation to 
obtain the resulting data word D R (k). If however, the strings 
represent packed data, the resulting data word D R (k) is, 
according to a preferred embodiment of the invention, derived by 
looking up a pre-generated value in a table 1310 based on the 

20 bit pattern given by the respective despread symbol string 
A, P (k), A| E (k), An_(k), AQp(k), A QE (k) and A QL (k). 

According to a preferred embodiment of the invention, in addition 
to the individual data words D R . )E (k), D R .| P (k), D R . )L (k), D R . QE (k), 
D R . QP (k) and D R . QL (k), corresponding accumulated data words are 
25 produced, such that after having..generated a data word D R (k), the 
sum of all data words D R (1) to D R (k) is also obtained. Hence, 
when the data word d(N) has been processed, the resulting sum 
D R (1) to D R (N) is also at hand. 

Moreover, the resulting data words represent payload information. 
30 For instance, based on the set of resulting data words D R .| E (k), 
D R .| P (k), D R ., L (k), D R . QE (k), D R . QP (k) and D R . QL (k) a demodulated 
piece of information of a transmitted data symbol sequence D as 
shown in figure 4 may be derived. In a steady state operation, 
basically only one of the resulting data words, say D R .| P (k), 
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actually carries the information. 

According to a preferred embodiment of the invention, the multi- 
plications involved in the multiplications between the interme- 
diate-frequency-reduced information words S| F .|(k) and S )F .Q(k) 
5 respectively and the modified code vectors CV m . P (k), CV m . E (k); 
CV m . L (k) are all performed by means of XOR- or SIMD-ope- 
rations. This is possible if the intermediate frequency infor- 
mation words SiF-i(k); S| F . Q (k) and the modified code vectors 
CV m . E (k), CV m . P (k); CV m . L (k) have compatible data formats. 

10 It was mentioned above with reference to figure 11 that more 
than one pair of early- and late pointers may be assigned to the 
prompt pointer in order to accomplish an improved tracking. 
Figure 14 illustrates how two pairs of such pointers P E1 ; P u and 
Pe2I Pl2 respectively are assigned around a prompt pointer P P 

15 according to an embodiment of the invention. The horizontal 
axis here represents a code shift CS and the vertical axis shows 
a normalised correlation factor. An optimal positioning of the 
prompt pointer P P is thus equivalent to a normalised correlation 
value of 1. The more the sampled data is shifted in either 

20 direction from this position, the lower the correlation value 
becomes down to substantially zero for a shift of one chip or 
more. The algorithm is perfectly balanced when (as illustrated in 
the figure) both the pointers P E i and P L1 of a first pair produce a 
first shifted correlation result, both the pointers P E2 and P U2 of a 

25 second pair produce a second shifted correlation result, the first 
shifted correlation result is larger than the second shifted corre- 
lation result and both the first and the second shifted correlation 
results are smaller than the correlation result obtained at the 
prompt position P P . Said correlation results are preferably 

30 generated according to a process corresponding to what has 
been described above with reference to figure 13. Based on the 
correlation results for the pointer positions P P , P E1 , P L i, P E 2 and 
P L2 a (if necessary new) set of modified code vectors CV m . E (k), 
CV m . P (k); CV m . L (k) is selected, which is deemed to be optimal. 
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Figure 15 shows a signal receiver 1500 for receiving navigation 
data signals transmitted in a navigation satellite system accor- 
ding to an embodiment of the invention. The receiver 1500 
includes radio front end unit 1510, an interface unit 1520 and a 
5 digital processor unit 1530. 

The radio front end unit 1510 is adapted to receive a continuous 
radio signal S H f. and in response thereto produce a corre- 
sponding electrical signal S| F which has a comparatively high 
frequency. The interface unit 1520 is adapted to receive the 

10 electrical signal S| F , and in response thereto, produce a 
sequence of sample values that represents the same information 
as the electrical signal S| F and is divided into data words d(k). 
The digital processor unit 1530 is adapted to receive the data 
words d(k), and in response thereto, demodulate a data signal. 

15 The digital processor unit 1530, in turn, includes a memory 
means 1535, which is loaded with a computer program that is 
capable of controlling the steps of the proposed procedure when 
the program is run in the processor unit 1530. 

In order to sum up, the general method of processing spread 
20 spectrum signals according to the invention will now be 
described with reference to a flow diagram in figure 16. 

A preparation step 1600 pre-generates code vectors which each 
represents a signal source specific code sequence that is 
intended to be received (or at least possible to receive) and 
25 demodulated in the receiver. The step 1600 is performed before 
the signal reception is initiated. 

A step 1610 then receives a continuous signal of a compara- 
tively high frequency. A following step 1620 samples the 
continuous signal at a basic sampling rate, whereby a resulting 
30 sequence of time discrete signal samples is produced. Each 
sample is also quantised (either with a relatively low-resolution, 
such as with 1 bit per sample value, or with a relatively high 
resolution depending on the application and the number of data 
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bits per samples delivered from the radio front end unit 1510), 
such that a corresponding level-discrete sample value is 
obtained. Subsequently, a step 1630 forms data words from the 
sample values, where each data word includes one or more 
5 consecutive sample values. After that, a step 1640 carries out 
correlation operations between the information in the data words 
and at least one of the pre-generated representations of a signal 
source specific code sequence. The correlation step involves 
correlating at least each vector in a sub-group of the code 

10 vectors with at least one vector that has been derived from a 
data word. A following step 1650 produces data as a result of 
the correlation performed in step 1640. Then, a step 1660 
investigates whether the sampled sequence has ended, i.e. 
whether there are any more data words to process. If more 

15 sampled data to process is found, the procedure returns to the 
step 1640. Otherwise, the procedure loops back to the step 
1610 again for continued reception of the incoming signal. 
Naturally, such signal is preferably also received during 
execution of the steps 1620 to 1660. The sequential procedure 

20 illustrated in figure 16 is merely applicable to a specific received 
signal segment. Preferably, all the steps are actually 1610 to 
1660 are actually performed in parallel. 

The process steps, as well as any sub-sequence of steps, 
described with reference to the figure 16 above may be 

25 controlled by means of a programmed computer apparatus, such 
as a microprocessor located in a GNSS-receiver. Moreover, 
although the embodiments of the invention described above with 
reference to the drawings comprise computer apparatus and 
processes performed in computer apparatus, the invention thus 

30 also extends to computer programs, particularly computer 
programs on or in a carrier, adapted for putting the invention 
into practice. The program may be in the form of source code, 
object code, a code intermediate source and object code such 
as in partially compiled form, or in any other form suitable for 

35 use in the implementation of the process according to the 
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invention. The carrier may be any entity or device capable of 
carrying the program. For example, the carrier may comprise a 
storage medium, such as a ROM (Read Only Memory), for 
example a CD (Compact Disc) or a semiconductor ROM, or a 
5 magnetic recording medium, for example a floppy disc or hard 
disc. Further, the carrier may be a transmissible carrier such as 
an electrical or optical signal which may be conveyed via 
electrical or optical cable or by radio or by other means. When 
the program is embodied in a signal which may be conveyed 
10 directly by a cable or other device or means, the carrier may be 
constituted by such cable or device or means. Alternatively, the 
carrier may be an integrated circuit in which the program is 
embedded, the integrated circuit being adapted for performing, 
or for use in the performance of, the relevant processes. 

15 The term "comprises/comprising" when used in this specification 
is taken to specify the presence of stated features, integers, 
steps or components. However, the term does not preclude the 
presence or addition of one or more additional features, 
integers, steps or components or groups thereof. 

20 The invention is not restricted to the described embodiments in 
the figures, but may be varied freely within the scope of the 
claims. 
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Claims 

1. A method of processing spread spectrum signals com- 
prising: 

receiving a continuous signal (S| F ) of a comparatively high 
5 frequency, 

sampling the continuous signal (S| F ) at a basic sampling rate 
(r s ), whereby a resulting sequence of time discrete signal 
samples (S[si]) is produced, 

quantising each signal sample into a corresponding level- 
10 discrete sample value, 

forming of a plurality of data words (d(1) d(N)) which 

each includes one or more consecutive sample values (s 1( .... s n ), 
and 

correlation between information in the data words (d(k)) and 
15 at least one representation (CS(i)) of a signal source specific 
code sequence (CS), characterized by 

the method comprising a preparation for the correlation step, 
wherein, before receiving the continuous signal (S| F ), a multitude 
of code vectors (CV; CV m ) are pre-generated, each code vector 

20 representing a particular code sequence (CS(i)) of the at least 
one signal source specific code sequence (CS), and 

the correlation step involving multiplying at least each vector 
in a sub-group (CV m . E , CV m . P ; CV m . L ) of the code vectors (CV; 
CV m ) with at least one vector (S )F -i(k); S| F -o(k)) derived from the 

25 data word (d(k)). 

2. A method according to claim 1, characterized by each 
code vector (CV; CV m ) representing a particular signal source 
specific code sequence (CS(i)) being sampled at the basic sam- 
pling rate (r s ) and quantised with the quantising process being 

30 used to produce the level-discrete sample val ues (si, s n ). 

3. A method according to any one of the preceding claims, 
characterized by the at least one signal source specific code 
sequence (CS) being pseudo random noise. 
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4. A method according to any one of the preceding claims, 
characterized by the receiving step involving down conversion 
of an incoming high-frequency signal (S HF ) to an intermediate 
frequency signal (S| F ), the high-frequency signal (S H f) having a 

5 spectrum which is symmetric around a first frequency (f H F) and 
the intermediate frequency signal (S )F ) having a spectrum which 
is symmetric around a second frequency (f| F ) being considerably 
lower than the first frequency (f HF ). 

5. A method according to claim 4, characterized by the steps 
10 of: determining a maximum frequency variation (f D ) of the 

second frequency (f| F ) due to Doppler-effects, 

defining a Doppler frequency interval (fiF-mm - fiF-max) around 
the second frequency (f| F ), the Doppler frequency interval (fiF-min - 
fiF-max) having a lowest frequency limit (fiF-min) equal to the 
15 difference between the second frequency (f| F ) and the maximum 
frequency variation (f D ), and a highest frequency limit (f| F . max ) equal 
to the sum of the second frequency (f !F ) and the maximum 
frequency variation (f D ), 

dividing the Doppler frequency interval (f| F . m j n - fi F -max) into an 
20 integer number of equidistant (Af) frequency steps (f| F . min , f| F . m j n + Af, 
fi F -max)i and 

defining a frequency candidate vector (f| F -c) for each fre- 
quency Step (f| F -min. f| F -min + Af f| F -max)- 

6. A method according to claim 5, characterized by the steps 
25 of: determining an integer number of initial phase positions 

(<P 0 ^7) for the frequency candidate vector (f| F -c), and 

defining a carrier frequency-phase candidate vector 
(V f(p (f| F -c, ^c)) for each combination of carrier frequency candi- 
date vector (f| F -c) and initial phase position (<P 0 ^7). 

30 7. A method according to claim 6, characterized by the steps 
of: determining the number of elements (s n ) in each carrier 
frequency-phase candidate vector (V f<p (f| F . c , <Pc)). and 
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storing the carrier frequency-phase candidate vector 
(V f( p(f|F-c. ^c)) according to a data format being adapted to per- 
forming multiplication operations between the data word (d(k)) 
and a segment of a carrier frequency-phase candidate vector 

5 (V f(p (f IF . c , 9c)). 

8. A method according to claim 7, characterized by the 
adapting of the data format involving adding at least one ele- 
ment to each segment of the carrier frequency-phase candidate 
vector (V f(p (fiF-c. 9 C )) such that the segment attains a number of 
10 elements which is equal to the number of elements in the data 
word (d(k)), thus enabling a segment of the carrier frequency- 
phase candidate vector (V f< p(f| F -c. ^c)) and one of the at least 
one vector (S| F -i(k); S| F .Q(k)) to be processed jointly by at least 
one of a SIMD-operation and an XOR-operation. 

15 9. A method according to any one of the preceding claims, 
characterized by the steps of: 

determining a maximum variation a code rate due to 
Doppler-effects, 

defining a Doppler rate interval (CR C -min - CR C -max) around 
20 a center code rate (CR C ), the Doppler frequency interval having 
a lowest, code rate limit (CR C - m in) equal to the difference 
between the center code rate (CR C ) and the maximum code rate 
variation (CR D ), and a highest frequency limit (CR c . ma x) equal to 
the sum of the center code rate (CR C ) and the maximum code 
25 rate variation (CR D ), 

dividing the Doppler rate interval (CR C -min - CR C -max) into 
an integer number of equidistant (ACR) code rate steps (CR C -min. 
CRc-min+ACR CRc-max) and 

defining a code rate candidate (CR C -c) for each code rate 

30 Step (CRc-mim CRc-min + ACR, CRc-max)- 

10. A method according to claim 9, characterized by deter- 
mining an integer number of possible initial code phase positions 
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(0.0 0.9) for each code rate candidate (CR C -c)- 

11. A method according to claim 10, characterized by 
defining, for each signal source specific code sequence (CS(i)), 
a set of combinations between code rate candidate (CR C -c) and 

5 code phase position (0.0 0.9), each combination repre- 
senting a code rate-phase candidate vector. 

12. A method according to claim 11, characterized by 
generating, for each signal source specific code sequence 
(CS(i)), a set of code vectors (CV) by: 

10 sampling each code rate-phase candidate vector with the 

basic sampling rate (r s ) whereby a corresponding code vector 
(CV) is produced. 

13. A method according to any one of the preceding claims, 
characterized by generating a modified code vector (CV m ) on 

15 basis of each code vector (CV) by: 

copying a particular number of elements (E e ) from the end 
of an original code vector (CV) to the beginning of the modified 
code vector (CV m ), and 

copying the particular number of elements (E b ) from the 
20 beginning of the original code vector (CV) to the end of the code 
vector (CV m ). 

14. A method according to claim 13, characterized by storing, 
for each signal source specific code sequence (CS(i)), a set of 
modified code vectors ({CV m (CR C -c. Cph)}), where 

25 each modified code vector (CV m ) contains a number of 

elements (St, .... s m ) representing a sampled version of at least 
one full code sequence (CS), 

a particular modified code vector (CV m ) is defined for each 
combination of code rate candidate (CRc-c) and code phase 

30 position (Cph). 
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15. A method according to claim 14, characterized by 
adapting the data format of the modified code vectors (CV m ) with 
respect to the data format of the at least one vector (S| F -i(k); 
SiF-o(k)) derived from the data word (d(k)), such that a modified 

5 code vector (CV m ) and one of the at least one vector (S| P .|(k); 
SiF-a(k)) may be processed jointly by at least one of a SIMD- 
operation and an XOR-operation. 

16. A method according to any one of the claims 14 or 15, 
characterized by involving an initial acquisition phase and a 

10 subsequent tracking phase, wherein during the acquisition phase 
a set of preliminary parameters are established which are re- 
quired for initiating a decoding of signals received during the 
tracking phase, a successful acquisition phase resulting in at 
least one signal source specific code sequence (CS) being identi- 
15 fied and a transmitted signal being rendered possible to track, 
each of the at least one signal source specific code sequence 
(CS) being associated with tracking characteristics in the form of: 
a modified code vector (CV m ), 
a carrier frequency candidate vector (fiF-c). 
20 an initial phase position (9c), 

a code phase position (Cph), and 

a code index (CI) denoting a starting sample value for the 
modified code vector (CV m ). 

17. A method according to claim 16, characterized by the 
25 tracking involving: 

calculating, based the tracking characteristics, a prompt 
pointer (P P ) for each modified code vector (CV m ), the prompt 
pointer (P P ) indicating a code sequence start position, and an 
initial prompt pointer (P P ) being equal to the code index (CI), and 
30 assigning, around each prompt pointer (P P ), at least one 

pair of early- and late pointers (P E , Pi_; Pel Pli Pe2. Pl2>. where 
the early pointer (P E ) specifies a sample value being positioned 
at least one element before the prompt pointer's (P P ) position, 
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and the late pointer (P L ) specifies a sample value being 
positioned at least one element after the prompt pointer's (P P ) 
position. 

18. A method according to claim 17, characterized by the 
5 tracking involving: 

receiving a sequence of incoming level-discrete sample 
values (1210), 

forming data words (d(1) d(N)) of the sample values 

(1210) such that each data word (d(k)) contains a number of 
10 elements which is equal to the number of elements (s n ) in each 
carrier frequency-phase candidate vector (V f(p (f|F-c. <Pc)). 

calculating a relevant set of carrier frequency-phase candi- 
date vectors (V f(p (fiF.c. ^c)) for the data word (d(k)), and 

acquiring, for each carrier frequency-phase candidate 
15 vector (Vf V (f| F -c. ^c)) in the relevant set, a pre-generated in- 
phase representation (f| F i) and a quadrature-phase represen- 
tation (f )FQ ) of the vector (V f<p (fiF.c. ^c)) respectively. 

19. A method according to claim 18, characterized by the 
tracking involving: 

20 multiplying each data word (d(k)) with the in-phase repre- 

sentation (f| F |) of the carrier frequency-phase candidate vector 
(Vf<p(fiF-c. ^c)) in the relevant set to produce a first intermediate- 
frequency-reduced information word (SiF-i(k)), 

multiplying each data word (d(k)) with the quadrature- 

25 phase representation (f ( Fo) of the carrier frequency-phase candi- 
date vector (V f<p (f|F- C . <Pc)) in the relevant set to produce a second 
intermediate-frequency-reduced information word (S| F .Q(k)). 

20. A method according to claim 19, characterized by the 
tracking involving: 

30 multiplying the first intermediate-frequency-reduced infor- 

mation word (SiF-i(k)) with a modified code vector (CV m . P (k)) 
starting at a position indicated by the prompt pointer (P P ) to 



WO 2004/036238 



38 



PCT/SE2003/001544 



produce a first prompt-despread symbol string (A| P (k)), 

multiplying the first intermediate-frequency-reduced infor- 
mation word (S )F .|(k)) with a modified code vector (CV m . E (k)) 
starting at a position indicated by an early pointer (P E ) to 
5 produce a first early-despread symbol string (A| E (k)), 

multiplying the first intermediate-frequency-reduced infor- 
mation word (SiF-i(k)) with a modified code vector (CV m . L (k)) 
starting at a position indicated by a late pointer (P L ) to produce 
a first late-despread symbol string (A| L (k)), 
10 multiplying the second intermediate-frequency-reduced in- 

formation word (S !F . Q (k)) with a modified code vector (CV m . P (k)) 
starting at a position indicated by the prompt pointer (P P ) to 
produce a second prompt-despread symbol string (A QP (k)), 

multiplying the second intermediate-frequency-reduced in- 
15 formation word (S| F .Q(k)) with a modified code vector (CV m . E (k)) 
starting at a position indicated by the early pointer (P E ) to 
produce a second early-despread symbol string (A QE (k)), and 

multiplying the second intermediate-frequency-reduced in- 
formation word (S| F . Q (k)) with a modified code vector (CV m . L (k)) 
20 starting at a position indicated by the late pointer (P L ) to 
produce a second late-despread symbol string (A QU (k)). 

21. A method according to claim 20, characterized by the 
tracking involving deriving, for each despread symbol string, a 
resulting data word (D R .,p(k), D R ., E (k), D R ., L (k), D R . QP (k), D R . QE (k); 

25 D R . QL (k)). 

22. A method according to claim 21, characterized by deriving 
the resulting data words (D R .| P (k), D R ., E (k), D R .| L (k), D R . QP (k), D R . 
QE (k); D R . QL (k)) by looking up a respective pre-generated value in 
a table (1310). 
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23. A method according to any one of the claims 19 - 22, 
characterized by performing the multiplication between the data 
word (d(k)) and the in-phase representation (f| F i) of the carrier 
frequency-phase candidate vector (V f(p (fiF-c. ^c)) respective 

5 between the data word (d(k)) and the a quadrature-phase repre- 
sentation (fiFd) of the carrier frequency-phase candidate vector 
(V f(p (fiF-c, <Pc)) by means of at least one of a SIMD-operation and 
an XOR-operation 

24. A method according to any one of the claims 19 - 23, 
10 characterized by performing the multiplication between the inter- 
mediate-frequency-reduced information words (S| F .i(k); S| F . Q (k)) 
and the modified code vectors (CV m . P , CV m . E ; CV m . L ) by means of 
at least one of a SIMD-operation and an XOR-operation. 

25. A method according to any one of the claims 19 - 24, 
15 characterized by propagating, in connection with completing 

the processing of a current data word (d(k)) and initiating the 

processing of a subsequent data word (d(k+1 )): 

a pointer (p d ) indicating a first sample value of the 

subsequent data word (d(k+1 )), 
20 a group of parameters describing the relevant set of carrier 

frequency-phase candidate vectors (V f(p (f| F .c, <Pc)). 
the relevant set of code vectors (CV m ), and 
prompt-, early-, and late pointers (P P , P E , Pi)- 

26. A computer program directly loadable into the internal 
25 memory of a computer, comprising software for controlling the 

steps of any of the claims 1-25 when said program is run on 
the computer. 

27. A computer readable medium, having a program recorded 
thereon, where the program is to make a computer control the 

30 steps of any of the claims 1 - 25. 
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28. A signal receiver (1500) for receiving navigation data 
signals transmitted in a navigation satellite system comprising: 

a radio front end unit (1510) adapted to receive a conti- 
nuous radio signal (S H f) and in response thereto produce a 
5 corresponding electrical signal (S| F ) comparatively high 
frequency, 

an interface unit (1520) adapted to receive the electrical 
signal (S| F ) and in response thereto produce a sequence of 
sample values being divided into data words (d(k)), and 
10 a digital processor unit (1530) adapted to receive the data 

words (d(k)) and in response thereto demodulate a data signal, 
and including a memory means (1535), characterized in that 
a computer program according to claim 26 is loaded into the 
memory means (1535). 
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D 
D 


03/06/03 








IIC 

uo 


CC9001 A 
D003O 14 


R 
D 


14/10/03 








IIC 

uo 


9nni nn99ftQ 


A 

A 


31/05/01 








us 


2001002458 


A 


31/05/01 








us 


2001047243 


A 


29/11/01 








us 


2002015439 


A 


07/02/02 








us 


2002143467 


A 


UO/ XV/ MC 








us 


2002146065 


A 


10/10/02 








us 


2002150150 


A 


17/10/02 








wo 


9740398 


A 


30/10/97 








us 


5901171 


A 


04/05/99 








us 


6041280 


A 


21/03/00 








us 


6522682 


B 


18/02/03 








us 


2003123531 


A 


03/07/03 


US 5592173 


A 


07/01/97 


NONE 
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